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ABSTRACT 

We present results of fitting the 50-day time series of photometry of a Cen A taken 
by the WIRE satellite in 1999. Both power spectrum and autocovariance function 
(ACF) fitting techniques were used in an attempt to determine mode frequencies, 
rotational splittings, lifetimes and amplitudes of low-i? p-modes. In all, using both 
techniques, we managed to fit 18 modes (seven I = 0, eight I = 1 and three I = 
2) with frequencies determined to within 1-2 /iHz. These estimates are shown to 
. be 0.6 ± 0.3 /iHz lo wer, on average, than the frequencies determin ed from two other 

' more recent studies l|Bouchv fc Carriei1l2002t IBedding et al.ll2004lh which used data 

\Q \ gathered about 19 months after the WIRE observations. This could be indicative of 

an activity cycle, although due to the large uncertainty, more data would be needed 
to confirm this. 

Over a range of 1700 to 2650 /iHz we were also able to use the ACF fitting to 

determine an average lifetime of 3.9 ± 1.4 days, and an average rotational splitting of 

0.54 ± 0.22 /iHz, which is the first ever reliable estimate of this parameter. In contrast 

to the ACF, the power spectrum fitting was shown to return significantly biased results 

52 ' for these parameters. 
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1 INTRODUCTION der to initially detect as many modes as possible and then 

to better constrain the limits placed on the determined fre- 
The past ten years nave seen a number of increasingly sue- TJ . r , . , . , , . 

, , , quencies. It is of course also important to improve resolution, 

cessful attempts to detect and measure solar-like oscillations 

in other stars. Due to its proximity and similarity to the Sun, 

many of these studies have been focused on the star a Cen A. 



but practical constraints have meant all these studies were 

limited in the length of observations that could be made. 

This has meant accurate determination of mode parameters 
The first clear detection of p-modo oscillations on this star , . . , ..... , ,. , , 

. . it—; ' r j ITT— r _ . such as power, rotational splitting and lifetime has been dit- 

was made bv Ibchou fc Bu zasi (2001) trom photometry us- g cu lt 

ing the Wide-Field Inf rared Explorer (WIRE) satellite taken 

over a 50-day period. ISchou fc Buzlsl J2OOIT1 correctly de- Here > we a PP ! y two sophisticated fitting procedures to 

termined the large frequency separation but, unfortunately, the WIRE 01 Cen A data collected in 1999 in order to im- 

wrong I identifications were made and hence an incorrect P rove the Parameter determinations. Although this data set 

value for the small separation was determined. Further de- has the Poorest SNR of the three studies mentioned above, 

tec tions and the first correct mode identifications were made 11 does have the longest time series. Hence, we would expect 

bv lBouchv fc Carrierl J2O0I) using a 13-day run of velocity to extract more reliable estimates of the average lifetime and 

measure ments taken by the C ORALIE spectrograph. More rotational splitting of the a Cen A modes. 



recently, IBedding et alJ ()2004Fl determined the frequencies The first fitting procedure we applied was a tradi- 

for over 40 individual modes from observations by the UVES tional power spectrum fitting method. This involved tak- 
and UCLES instruments taken over a period of 5 nights. ing the Fourier transform of the time series and then fitting 

The main driving force behind each subsequent study a Lorentzian-like model to the various mode peaks in the 
has been to improve the signal-to-noise ratio (SNR) in or- resulting power spectrum. The second procedure used is a 

new technique based on fitting the autocovariance function 
(ACF) of the time series (i.e., the unnormalized autocorre- 
* E-mail:stfletch@bison. ph. bham.ac.uk lation function). Since the modes seen in the ACF are all 
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superimposed, one must first filter the time series in order 
to isolate the modes one is hoping to fit. The ACF can then 
be computed and a model based on an exponentially de- 
caying, periodic function fitted to the result. This technique 
was first i ntroduced as a possib le method for mode determi- 
nation by|GabrieLel_alJ ll9_9_£[ ) and d eveloped more fully by 
iFletcher. Chaplin fc Elsworthl (12004) in an attempt to bet- 
ter constrain mode parameters of long solar p-mode data 
sets. 

In the init i al in vestigation of the WIRE data by 
ISchou fe BuzasJ i200lf) only a handful of modes were identi- 
fied due to the poor SNR in the data set. However, a distinct 
advantage of revisiting this data comes in having a large 
number of robust mode identifications from the aforemen- 
tioned two later studies. This provides additional a-priori 
information that we can use as initial 'guess' values for our 
fitting procedures. 

In order to fully test our fitting of the WIRE data we 
also generated a set of artificial time series. These data were 
created specifically to mimic the WIRE data allowing us 
to explicitly determine the precision and robustness of our 
fitted parameters. We detail the creation of this simulated 
data in Section [5] In Section [3] we describe in detail the 
procedure involved in applying the two fitting techniques 
to the data. Finally, in Section |1J we go on to present and 
analyse the results of our fitting for the mode frequencies, 
amplitudes, lifetimes and rotational splitting parameters. 



2 DATA 

The WIRE satellite collected 50 days of photometry obser- 
vations on the star a Cen A between July 15 and Septem- 
ber 3, 1999. The orbit was 96 minutes, low Earth and 
Sun-synchronous, which unfortunately lead to rather signif- 
icant pointing constraints. Hence, observations could only 
be made for about 40 minutes out of each orbit. Addition- 
ally there were some days when no data were collected and 
some days where data could not be used due to process- 
ing problems. Hence, the resulting duty cycle for the WIRE 
time series was only 15%. The processing of th e WIRE data 
to pro duce the final time ser i es is d escribed in lBuzasi et alJ 
feOOd) and ISchou fc Buzasil (1200 J) . The data were binned 
to a 10s cadence resulting in a series of 432,000 points. The 
processed time series is shown in FigQ for both the full 50 
days and a zoomed in section with just the first 8 days in 
order to show the main gap structure. Being as the time 
series was made using photometry observations, the ampli- 
tudes are given in parts-per-million (ppm). The method for 
converting between ppm and velocity is described in Sec- 
tion ED 

A full set of simulated data made to resemble the WIRE 
time series was produced in order to test fitting on the 
real data. The Laplace transform solution of the equation 
of a forced, damped harmonic oscillator was used in or- 
der to generate individual mode components of the artifi- 
cial time serie s. The applica t ion of this model is described 
more fully in IChaolin et ail (|l99/t l. A set of low-£ modes 
covering the ranges ^ I ^ 2 and 1700 ^ v ^ 2650 /iHz 
was created in this way. It should be noted that the £ = 3 
modes are not included since they do not appear to be de- 
tectable in the WIRE photometry data, although for more 




(a) 50 days 



(b) 8 days 

Figure 1. The Wire a Cen A time series over (a) the full 50 days 
and (b) the first 8 days only. 



recent velocity observations they are l)Bedding et alJl2004h . 
For each £ 7^ 0, a rotationally induced splitting pattern was 
included, for which the visibilities, e, of each component were 
fixed using the following equ ations l|Qizon fc Solankll2003l : 
iToutain fc Gouttebrozelll993ft 



£l,0 = COS I, 



Ei.ii = h sin 2 L 



£2,0 = j (3 cos 2 i - 1) 
£2,±i = | sin 2 (2i), 



3 • 4 • 

£2, ±2 = s sin I, 



(1) 

(2) 
(3) 
(4) 
(5) 
79° 



and the inclination, i, was chose n 
iPourbaix. Neuforge-Verheecke fc Noelslll999l) . 

A database of frequency and powers, based on prelim- 
inary fits to the WIRE data, was used in order to fix the 
various characteristics of each mode. Where no mode s could 
be identified, frequencies presented in lBedding et al] (l2004h 
were used. 

Based on fits to the ACF of the WIRE data (see sec- 
tion ^. 21 we chose a constant rotational splitting of 0.5 ^Hz 
and a constant linewidth (full width half maximum) of 1.0 
piHz (equivalent to a lifetime or e-folding amplitude time of 
~3.7 days). Analysis of fits to these artificial data (see sec- 
tions 4.2 and 4.3) showed the parameters could be robustly 
determined. This gives confidence that the simulated spec- 
tra were a good match to, and a reliable cross check for, the 
real data. A background offset (white noise) was also added 
so as to match the WIRE data. In total 50 independent data 
sets were created, each of a length equal to the WIRE time 
series. The WIRE window function was imposed on each of 
them. 
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Table 1. n and t values included in fitting. 
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(a) The WIRE Power Spectrum 



Frequency (pHz. 



(b) The WIRE Spectral Window 

Figure 2. A (a) The WIRE Power Spectrum over the region 1500 
SC v sC 2850 /iHz and (b) the spectral window plotted over the 
same total frequency range. 



3 MODELLING AND FITTING 

In this section we outline the two different methods used to 
determine the various p mode parameters. First we discuss 
the more traditional power spectrum technique and then go 
onto introduce the ACF method. In both cases we explain 
the specific models used to fit the WIRE data. 



3.1 Power Spectrum 

When fitting the power spectrum of a solar or solar-like p- 
mode time series the various mode components are generally 
modelled using either a Lorentzian, or an asymmetric func- 
tion that models small depa rtures from a symmetric shape 
(Niga rrTfe Kosovichevlll99M) . For low-£ (Sun- as- a- star) data 
this is commonly done by fitting individual mode pairs (i.e., 
£ — with neighboring £ — 2; and I = 1 with £ — 3) using 
narrow fitting 'windows' centered on the target pair. 

Unfortunately, certain characteristics of the WIRE data 
make this strategy difficult to implement. The rather severe 
window function mentioned in Section[3]results in prominent 
sidebands at 173.6 ^iHz from the mode peaks. As a result 
the sidebands of modes lie in the vicinity of other overtones. 
This makes it very difficult to isolate a single mode pair and 
its sidebands without the presence of intervening modes. 
Therefore, we adopt a strategy of fitting Lorentzian profiles 
to all modes simultaneously. Note that the WIRE power 
spectrum and the spectral window are shown in Fig. |5] 

In the 1700 ^ v ^ 2650 /iHz range a total of 27 modes 
were included in the simulated data sets: 9 each of £ = 0, 1 
and 2. The £ — 1 and 2 modes were split into three and five 



components respectively to match the rotational splitting 
pattern. However, we chose to fit only the modes that are 
identifiable in the WIRE data, of which there are eighteen. 

For the power spectrum method, the power spectral 
density, P, was modelled using a standard Lorentzian profile 
for each peak summed over all visible modes and correspond- 
ing sidebands in the chosen frequency range, offset with a 
flat background, B, i.e., 



P — B + H n lmkZ n l mk , 



(6) 



where 



v — v n i + ras + kw 
A/2 



Here, H is the height of each mode component (i.e., max- 
imum power spectral density), v n i is the central frequency 
of a mode, s is the rotational splitting, w is the sideband 
spacing and A is the width. 

In the WIRE data, over the range 1700 < v < 2650 
/xHz, the radial overtone number, n, increases from 14 to 
23. However, since we fit only 18 modes we do not sum over 
all the possible n£ combinations. Table 1 gives the list of 
n and ^-valued modes that are fitted, m is the azimuthal 
order and is summed over — £ ^ m ^ £, while k is a dummy 
variable allowing us to include the first order sidebands in 
our model and as such is summed over — 1 ^ k ^ 1. 

Sidebands were assumed to lie at a fixed spacing of 173.6 
piHz and to have the same width as the main peaks. The 
fractional height of the sidebands were also fixed according 
to the ratio of the sidebands to the main peak in the fourier 
transform of the window function. In this way, we did not 
fit any parameters associated solely with the sidebands. 

The number of fitted parameters was further reduced by 
assuming fixed height ratios for the rotationally split compo- 
nents. These were fixed according to equations Q to |S] again 
assuming an inclination for a Cen A of 79°. 

Finally, it should be noted that because of the relatively 
short duration of the time series and the poor SNR, we did 
not feel that the data justified the use of the more compli- 
cated asymmetric iNigam fc Kosovichevl (I1998T) model. 

To fit the model to the power spectrum we used a 
Powell multi-dimensional hill-climbing minimization algo- 
rithm, which maximized an appropriate log-likelihood func- 
tion. This function was based on the assumption that the 
power spectrum is distributed with negative exponential 
(i.e., x 2 , with two degrees of freedom) statistics. The follow- 
ing parameters were varied iteratively until they converged 
on their best fitting values. 

(i) A frequency for each mode. 
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(ii) A single height for each mode. Heights of rotationally 
split components were fixed relative to the strongest outer 
I = \m\ components (see equations Q to 0. 

(iii) A single width for all modes. 

(iv) A single splitting for the 1=1 and 2 modes 

(v) A flat, background offset for the whole fit. 

A total of 2M+3 parameters were therefore fitted where M 
is the number of modes included in the fitting (i.e, f 8 in this 
case). We should add that in order to recover the power, 
(i.e., the square of the amplitude), A 2 , in each mode from 
the fitted height, H, and the width, A, we used the following 
expression: 



A 2 = H (|TA+ 1 



(7) 



where T is the length of the time series. Details on the 
derivation of this expression are given in the appendix. 



3.2 The Auto-Covariance Function 

The Autocovariance Function (ACF) is the product of a 
data series with a shifted versions of itself over successive 
time lags, r. For a time series of discrete measurements, Xi, 
it is defined as: 



E< 



(JV-l)- 




(3<i CC^i^CCi 



x) 



N - 



(8) 



where TV is the number of points in the data series and x 
is the mean of the sample. In many data series the mean is 
often very close to zero, as is the case for p-mode intensity or 
Doppler velocity residuals, in which case the ACF reduces 
to: 



V-^(.JV- 



(JV-l)- 



N 



(9) 



The more commonly used normalized version of the au- 
tocovariance function is termed the autocorrelation function 
and for a time series with zero mean is given by: 



Pt = 



El 



(iV-1)- 




XiXi-\-T 



N - 



•EI 



(io) 



However, by fitting the autocovariance function we can di- 
rectly obtain absolute estimates of the powers associated 
with the modes, as opposed to relative powers that would 
be obtained from fitting the autocorrelation function. 

Since the ACF is computed in the time domain, the 
periodic waveforms of the various components are superim- 
posed. Therefore, to fit a certain set of modes in a given 
frequency range one must first apply a band-pass filter to 
the time series. As the observable modes in the WIRE data 
lie in the region 1700 ^ v 2650 /iHz we need to at least 
filter over this range. However, we must also take into ac- 
count the effect of the window function. As one would ex- 
pect, the sidebands in the power spectrum manifest in the 
ACF as additional waveforms. Hence, we choose to extend 
the bandpass range to 1500 ^ v sC 2850 /iHz in order to 
include all the sideband frequencies. 

In Fig.[3]we show the normalized ACF (i.e., the autocor- 
relation function) of the WIRE time series, filtered between 
1500 ^ v ^ 2850 /^Hz, and plotted over: the first 3600 points 
(10 hours) in Fig. 13(a)] and over the first 43200 points (120 



(a) ACF - 10 hours 




(b) ACF - 120 hours 

Figure 3. Normalized ACF of WIRE data over (a) the first 10 
hours and (b) the first 120 hours. The scale has been reduced in 
(b) in order to better show the features. 



hours) in Fig. |3(b)| Th e filtering was p erformed using a non- 
recursive digital filter ljWalravenlll984h . There are a number 
of distinctive features in the ACF the most striking being 
the very sharp decay in the ACF that occurs over the first 
few hours. This is due to the poor SNR in the WIRE data 
and is explained in more detail below. 

After this initial decay the ACF structure is dominated 
by the p-mode signal. On the smallest timescale, there is a 
quasi sinusoidal waveform with a period of about 7.5 min- 
utes. This particular feature can be seen in Fig. |3(a)| but 
is too rapid to be seen in Fig. |3(b)| It is due to the main 
set of modes that are seen at around 2000 /iHz in the cv-Cen 
spectrum. On a longer timescale we see quasi-periodic struc- 
ture close to ~1.5 hours. This is due to the prominent side- 
bands which occur at 173.6 /iHz to either side of the main 
peaks. We also see a hint of a ~5-hour quasi-periodic struc- 
ture due to the separation between 1 = and £ = 1 modes. 
Longer time scale quasi-periodic variations due to smaller 
mode spacings, such as those between adjacent modes in 
the low-£ pairs, and rotational splitting, are more difficult 
to observe and therefore extract. However, there is evidence 
of an overall exponential decay in the structure of the ACF 
due to the finite lifetimes of the modes. 

The ACF was modelled using a damped harmonic os- 
cillator equation summed over all visible modes and corre- 
sponding sidebands in the chosen frequency range with an 
additional background function, B T , 



f T = B T + Anlmk COS (iJnlmkti-) exp(—jt T ), 



(11) 



where, 
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ik = \J (2-7T (i/„i +ms + kw)) 2 



A 2 is the mode power, uj is the natural angular frequency 
of the mode, 7 is the damping constant and t is time. Note 
that 7 is related to the width of mode peaks by A = 7/71". 
B r is the background component but in the ACF this is 
not a simple offset but must be treated according to the 
expression: 



SUl(2lTV\tr) 
2lTVltr 



COs{2-KV2tr) 



where b is the power in the background and is the parameter 
to be fitted. v\ is a value given by the extent of the filtered 
frequency range divided by 2 (i.e., 675 //Hz in this instance) 
and V2 is the central frequency in the filtered range (i.e., 
2175 fiRz). The sine term is a direct result of filtering over a 
finite frequency range. The wider this range the higher the 
frequency and the faster the background function decays. It 
is this decay that dominates the first few time samples of the 
WIRE ACF. The extra waveforms due to the window func- 
tion are included in the model in the same way as sidebands 
were treated in the power spectrum. 

For the ACF, a gradient-expansion algorithm was used 
to perform a non-linear least-squares fit to our model, with 



X 



E 



[Yr-f(U, 



(12) 



the quantity to be minimized, where Y are the data, / is 
the model, t is time, a is the vector of parameters and a is 
the expected error on each point of the ACF. Note that we 
are essentially fitting the same set of parameters as with the 
power spectrum except the power (square of the amplitude) 
in the modes is fitted directly. 

We note that two assumptions are made in order to sim- 
plify fitting the model to the ACF. Firstly, we take the error 
distribution, a T to be constant over the range of the ACF 
that we fit. A plot of the standard deviation over all the 
ACF's of the simulated data shows this to be a reasonable 
approximation. Secondly, we ignore the effect of correlation 
between one point in the ACF and the next. While correla- 
tion clearly must be present, the effect of ignoring it can be 
shown to only affect attemp ted error ca lculations and not 
the fitted values themselves <lKuanll2000l) : hence the reason 
for using Monte Carlo simulations to fix the errors on our 
fits. 



4 RESULTS AND ANALYSIS 

In this section we present the results and analysis of fitting 
the WIRE data using the two methods described in Sec- 
tion [3] We initially show a graphical representation of the 
fitted model overlaid on the power spectrum before analyz- 
ing the four different types of parameters fitted (frequen- 
cies, linewidth, amplitudes and background) in more detail 
in separate subsections. 



4.1 Graphical Representation of the Fits 

In order to give a clear picture of the fits as determined from 
the power spectrum fitting technique, Fig.0]shows the fitted 



model overlaid on the WIRE power spectrum in the region 
1700 ^ 1/ ^ 2500 /iHz. This plot shows a very good illustra- 
tion of the fitted splittings. (Note that in sections 14.31 and 
I4.4l we show that fits to the power spectrum have a tendency 
to slightly overestimate the splitting and underestimate the 
width.) 

The fitted peaks that are not labelled in the plot are the 
first order sidebands. The dotted line extending outwards 
from the strong n=19, l=\ peak shows the frequency spac- 
ing at which the associated sidebands are placed. Having 
sidebands occur at such a large distances from their main 
peak is somewhat unusual for those more used to dealing 
with diurnal sidebands. 

On closer inspection of the plot, there are some fitted 
modes that have very small powers (namely the n=17, £=1, 
n=20, £—0 and the n=21, £=1 modes) and as such, the re- 
liably of these fits is somewhat questionable. The reason for 
their presence is that the fitting routine attempts to fit peaks 
for all the modes included in the model. Hence, if a mode 
is expected at a certain frequency, and there is no large and 
obvious peak in the nearby vicinity, the fitting routine may 
latch on to a smaller peak which is only associated with the 
background noise. However, with these fitted modes having 
such small powers, even if they are not true signatures of 
mode excitation, their effect on the overall width and rota- 
tional splitting estimates is negligible. 

In Fig.|S] four 15 piHz slices of the spectrum are shown, 
each located on fairly strong modes. For the £=\ and £=2 
modes, it is possible to identify a rotational splitting pattern 
in the spectrum, and the model is seen to fit reasonably well. 
While we would also like to be able to show similar graphical 
illustrations of the splitting and widths as determined from 
the ACF fitting, the superimposed nature of the waveforms 
makes this difficult. 



4.2 Mode Frequencies 

Both methods produced well-constrained frequency esti- 
mates for seventeen out of eighteen of the modes we tried to 
fit. Sixteen of these were determined via both methods while 
the other two modes were fitted by one method only. The 
fitted frequencies are shown in Tables |5] and |3] respectively 
where the quoted uncertainties were determined from the 
standard deviation of fits made to the 50 simulated spectra. 
The errors are quite large for some of the estimates due to 
the high background and resulting poor SNR of the WIRE 
data. A comparison between the two sets of results shows 
the differences to be very small. The root-mean-square (rms) 
difference is 0.3 //Hz, compared with an average error of 1.4 
/iHz. 

We can also compare our fre quenc ies with those de- 
termin ed by iBedding et alj (120041) and iBouchv fe Carrier! 
( 2002). Looking at individual modes we see differences of the 
order of 1 /iHz which is within the errors given. However, 
when taking a weighted average of the differences we find 
the estimates from the WIRE data are about 0.8 ± 0.3 £iHz 
lower than the frequencies determin ed from IBedding et alJ 
l|2004l) and IBouchv fc Carrier! J2002I) . One contribution to 
this difference comes from the shift in the mode frequencies 
due to the relative motion of the Earth around its orbit. By 
considering the times of year during which the observations 
were carried out, and the ecliptic coordinates of a Cen A, it 
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Figure 4. WIRE Spectrum for the region 1700 $C v ^ 2500 /xHz, with the fitted model overlaid . The fitted modes arc labelled and the 
horizontal lines extending from the strong peak at ~ 2200 //Hz indicate the scale of the sideband distance. Three of the modes have very 
questionable fits, those at n=17, i=l, n=20, £=0 and n=21, i=l. 
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Figure 5. 15 /iHz slices of the WIRE spectrum showing the fits for four different modes. 
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Table 2. Frequencies extracted by power spectrum fitting 
method (/iHz). Brackets indicate modes with questionable fits, 
see Section IQ 



n 


e=o 


1=1 


1=2 


14 






1725.4 ± 1.5 


15 




1775.2 ± 1.3 




16 


1838.2 ± 1.2 


[1886.6 ± 1.7] 




IV 


1942.4 ± 1.0 


1992.7 ± 1.8 


2037.9 ± 1.4 


18 




2096.4 ± 1.6 




19 


2152.5 ± 1.0 


2201.8 ± 0.7 




20 


[2257.2 ± 1.4] 


2307.1 ± 1.2 




21 


2363.2 ± 1.0 


[2414.6 ± 1.4] 


2465.5 ± 1.8 


22 


2469.0 ± 0.9 






23 




2623.9 ± 1.3 





Table 3. Frequencies extracted by ACF fitting method (/tHz). 



1600 1800 2000 2200 2400 2600 2800 
Frequency (/iHz) 



(a) Power Spectrum 



1600 1800 2000 2200 2400 2600 2800 
Frequency (fiHz) 



(b) ACF 



Figure 6. Results of fitting artificial WIRE-like data. Aver- 
age difference between fitted and input frequencies for £=0 (di- 
amonds), 1=1 (triangles) and 1=2 (squares) modes. Error bars 
indicate the error on the mean as given by the standard devia- 
tion in the fitted frequencies divided by the square root of the 
number of fits. 



n £=0 1=1 1=2 



14 






1725.3 ± 1.8 


15 




1774.9 ± 1.5 




16 


1838.5 ± 1.2 


1886.6 ± 1.7 




17 


1942.4 ± 1.3 


[1993.1 ± 1.8] 


2038.8 ± 1.9 


18 




[2096.9 ± 1.9] 




19 


2152.5 ± 0.9 


2202.1 ± 0.7 




20 


[2257.1 ± 1.5] 


2307.7 ± 1.1 




21 


2363.4 ± 0.8 




2465.6 ± 2.2 


22 


2469.2 ± 0.8 






23 


2578.3 ± 1.7 


2624.3 ± 1.4 





can be shown that this effect should lead to the WIRE data 
having frequencies ~ 0.2 /iHz lower than the other two data 
sets. Therefore, the corrected difference is 0.6 ± 0.3 /iHz, 
which is only significant at 2 sigma. 

Since the WIRE data set is taken 19 months prior to 
the others one might conjecture that the difference is due 
to an activity cycle on a Cen A. Were a 0.6 /iHz shift over 
a period of just 19 months to be real, it would suggest a 
large and rapid (or possibly just very large) activity cycle. 
By comparison, the minimum to maximum change in mode 
frequencies due to the solar cycle, is ~0.4 /iHz and this oc- 
curs over a 5.5-year period. Of course, a zero or very small 
shift is not ruled out and so it is clear that further and bet- 
ter quality sets of observations would be needed in order to 
make any solid conclusions. 

As well as using the simulated data to estimate the 
uncertainties on our fits we also used them to test accuracy 
and reliability. The plots in Fig. |5] show the mean difference 
between the fitted and input frequencies of the simulations. 
The associated error bars are the error on the mean, given 
by the standard deviation divided by the square root of the 
number of fits. For both methods the estimates are generally 
accurate to within errors and there does not seem to be any 
particular bias towards an under or overestimation of the 
input frequency. 

The errors on the real data estimates given in Tables H 
and|3]and on the artificial data estimates in Fig.HJshow how 
the uncertainties vary quite substantially from one mode 
to the next. This is to be expected since the uncertainties 



depend st rongly on the SNR, when the SNR is poor. T he 
equation jLibbrechtJll992l ; Flbutain fc Appourchauxlll994l) : 

<« = -T^V^+U^+l + V^) 3 > (13) 

gives an estimate of the error on the fitted frequency, o Vnl 
as a function of the linewidth, A, the length of the time 
series, T and the noise-to-signal ratio, j3 (i.e the inverse of 
the SNR). When the SNR is high, j3 is small, and the level of 
uncertainty will be dominated by A and T and as such will 
remain fairly constant from mode to mode. However, for a 
poor SNR, the p term is important and so small changes in 
this, due to the difference in amplitudes of the modes, will 
have a large effect on the uncertainty. 

The relationship between the SNR and the errors on the 
frequencies can be seen more clearly by plotting the stan- 
dard deviation against the input amplitude of the modes, 
as in Fig. Q (this is equivalent to plotting against the 
signal-to-noise in amplitude since the background is kept 
constant with frequency). The large SNR in the WIRE 
data also means that the typical uncertainties on our fre- 
quency estima t es are somewhat larger than those given by 
iBedding et"aH <2004ft . which were ~0.3 /iHz. 

4.3 Rotational Splitting 

There have to date been no successful attempts to mea- 
sure the rotational splitting of the p-modes in a Cen A and 
thereby make a first step to constraining the internal ro- 
tation. However, by applying the ACF fitting technique to 
the WIRE data and assuming a constant splitting across 
the frequency range we have been able to determine a rea- 
sonably well constrained estimate of 0.54 ± 0.22 /iHz. The 
simulated data were again used in order to check the accu- 
racy of this estimate. By fitting the ACF of the 50 artificial 
data sets an average splitting of 0.44 ± 0.03 /iHz was re- 
turned. This was within 2 standard errors of the input value 
of 0.5 /iHz, indicating that the fits are relatively unbiased. 
It should be noted that of the 50 estimates returned via the 
ACF method, 10 gave a splitting very close to zero. These 
'null' splitting results are a common problem when fitting 
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Table 4. Amplitudes extracted by power spectrum fitting method 
(cms" 1 ). Brackets indicate modes with questionable fits. 
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Figure 7. Results of fitting artificial WIRE-like data. Standard 
deviation of fitted frequencies as a function of the input amplitude 
of the modes. Symbols have the same meanings as in Fig. [H] 



poorly defined modes (either due to poor resolution, low 
SNR or a large width-to-splitting ratio), and we have taken 
the usual step of disregarding these values when determining 
the m ean and standard de viation of the fits to the simulated 
data iChaplin et al.ll200lTl . 

A value of the splitting was also returned when fitting 
the power spectrum. This method gave an estimate of 0.64 
± 0.25 /iHz, which is within errors of the estimate given by 
fitting the ACF. However, fits to the simulated data showed 
a tendency for the estimates of the splitting to overestimate 
the true input value. An average of 0.68 ± 0.03 was re- 
turned which gives a result which is ~6 standard errors too 
large. This is also a common problem in the solar case but is 
generally due to a large width-to-splitting ratio at high fre- 
quencies rather than a poor SNR. Because of the bias in the 
power spectrum estimate we would suggest that the value of 
the rotational splitting returned via fitting the ACF is more 
robust. 

A rotational splitting of 0.54 ± 0.22 /iHz gives a rota- 
tional period of 2llg 7 days. This agrees with a me asurement 
of the a Cen A surface rotation rate determined bv lJav et alJ 

(1996) of 23^2 days using chromospheric and transition re- 
gion emission line features as markers of ro tation. Also, a 
surfac e rotation rate of 2700 ± 700 ms _1 from lSaar fc Ostenl 

( 1997) coupled with theoretical predict ions of the alpha Cen 
A radius of 1.23 R {Morel et alJl200f^ . give a prediction of 
23 ± 6 days. A similar correspondence of the \cw-l splittings 
and the surface rate of rotation is seen for the Sun. 



4.4 Linewidths and Lifetime 

Since our model also assumed a fixed value for the linewidth 
we were able to obtain a well constrained estimate of this 
parameter. However we lose any knowledge of how the pa- 
rameter varies with frequency. Our estimate for the average 
linewidth of the modes from fitting the ACF of the WIRE 
time series is 0.92 ± 0.30 /iHz, which equates to a lifetime 
of 3.9 ± 1.4 days. Fits to the ACF of the simulated data 
gave an average width of 1.10 ± 0.05 /iHz, which is again 
within 2 standard errors of the input value of 1.0 /iHz. Our 
fits to the artificial data also showed a clear anti-correlation 
between the fitted linewidths and the fitted rotational split- 



n 


1=0 


i=\ 


1=2 


11 

15 
16 
17 


23+Jl 


13 +18 

[ 4±« ] 


41+34 
— 14 

23+_ 2 ° 


18 
19 
20 
21 
22 
23 


24±2 6 

[ 3±» ] 

38+ 1 * 

29+ 34 
-13 


13 +19 

76 +B4 
-30 
,, + 38 
d0 -14 

[ T±» ] 

26+" 
— 6 


24 +15 



tings. This meant that the fits that returned 'null' splittings 
also returned artificially high linewidths, as one might ex- 
pect. We therefore removed these fits from our calculation 
of the mean and standard deviation of the artificial results, 
in the same manner as for the rotational splitting. 

Our estimates of the linewidth from fitting the power 
spectrum gave us a much smaller value of 0.46 ± 0.38 /iHz, 
which is a lifetime of 8.1 ± 6.8 days. However, fits to the 
simulated data showed estimates of the linewidth to be sig- 
nificantly underestimated. An average fitted value across the 
50 data sets of 0.62 ± 0.06 /iHz was returned indicating 
an underestimate of over 6 standard errors. The same anti- 
correlation between linewidth and rotational splitting values 
that was seen when fitting the ACF was also seen here with 
overestimates of the splitting being combined with underes- 
timates in the linewidth. This leads us to conclude that the 
ACF fitting technique also gives a more reliable estimate of 
the linewidth and hence lifetime of the a Cen A p-modes. 

In relation to our assumption that fitting a constant 
width returns the average linewidth across the fitted modes 
we note that the ACF fitting has also been tested on artificial 
data simulating solar oscillations, for which the linewidth 
varies with frequency. In this case we found that fitting a 
constant width did indeed return the input average across 
the mod es. 

The lBouchv fc Carrier! J2002t and lBedding et alJ d2004l) 
studies were both made from observations lasting just a 
few days (5 and 12 respectively) and a direct measure- 
ment of the mode life time was difficult to achieve. However, 
iBedding et alJ <l2004l) were able to develop a method of es- 
timating the lifetime from the scatter in the frequency mea- 
surements. Their initial estimate put the average mode life- 
time at 1.4+0 4 days at 2.1 mHz, however, that has re cently 
been reevaluated as 2.3+J'g days llKieldsen et all20 05). This 
estimate is smaller than our value determined from fitting 
the ACF of the WIRE data, although the error bars overlap. 



4.5 Mode Amplitudes 

We have chosen to present the strength of the modes in 
terms of their amplitudes rather than power. We have also 
converted the units to velocity (cms -1 ) even though the 
WIRE observations were intensity measurements. In both 
cases this was done so as to more easily make comparisons 



Frequency, splitting, linewidth and amplitude estimates of low-l p modes of a Cen A 9 



Table 5. Amplitudes extracted by ACV fitting method (cms 1 ). 
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Figure 8. Results of fitting artificial WIRE-like data. Average ra- 
tios of fitted to input amplitudes. Symbols have the same meaning 
as m Fig. El 



with the iBedding et al.l i2004t) study. To convert between 
intensity given in parts pe r million (ppm) and velocit y we 
use the expression given in iKield sen fe Bedding! lll995f) : 



(¥). 



6L 



\550nm7 V5777K7 V 20.1ppm 



(14) 



osc denotes the amplitude of velocity oscillations in 
, (SL/L)x the change in intensity at the effective wave- 



where v, 
ms 

length A (~ 450nm for WIRE), T c g the effective temperature 
of the star (~ 5770 for a Cen A) and 5777K is taken as the 
effective temperature of the Sun. 

Of the four types of parameter examined in this paper 
the amplitudes are the least well constrained. This is be- 
cause we have attempted to fit a separate power to each 
mode, rather than an average value across several modes 
as was done with the linewidth and splitting. This results 
in large uncertainties on the amplitude estimates as shown 
in Tables [I] and |^| for the power spectrum and ACF fit- 
ting respectively. Indeed the errors on some of the weaker 
modes are larger than the actual estimated powers. Hence, 
we cannot put much credence in these values as they will be 
dominated by background noise. 

We have again used simulated data to test the accu- 
racy of the fitted powers. In a similar fashion to that used 
for the fitted frequencies, we have plotted in Fig. [S]the dif- 
ference between the power estimates averaged over the 50 
artificial data sets and the input values. The plots show 
that both methods return a number of estimates that sig- 
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Figure 9. Results of fitting artificial WIRE-like data. Average 
ratios of fitted to input amplitudes, plotted as a function of input 
amplitude. Symbols have the same meaning as in Fig. IB1 



nificantly underestimate the true input amplitude. In Fig. 
we plot the same data but as a function of the mode ampli- 
tude. This shows that when fitting the ACF, the amplitude 
of the stronger modes seem to be underestimated by around 
30 percent, whereas the weaker modes suffer less bias. This 
pattern is less obvious when fitting the power spectrum with 
little evidence of the bias being dependent on the strength 
of the modes. Fig.|H| also shows a tendency for the power of 
the lower ^-valued modes to underestimate the input values 
more so than for higher ^-valued modes. This is particulary 
noticeable for the power spectrum fits. 



5 SUMMARY 

The 50-day time series of photometry observations taken 
in 1999 by the WIRE spacecraft has been reanalysed using 
power spectrum and autocovariance fitting methods. With 
the help of a-priori information, regarding the location in 
frequency of modes from other recent a Cen A studies, 
we have managed to fit 18 different modes in the power 
spectrum and autocovariance function (ACF), 16 of which 
were fitted in both. The values of the fi tted frequencies are 
slig htly lower than those det ermined bv lBedding et afll2004l 
and lBouchv fc Carrier! l200l although without better data 
we cannot say with any confidence whether this is indicative 
of an activity cycle for a Cen A. 

In addition to the frequencies we have also been able 
to estimate an average rotational splitting across the fitted 
modes of 0.54 ± 0.22 ^tHz using the ACF. An average life- 
time was also estimated by fitting the ACF and was found 
to be 3.9 ± 1.4 days. Although th e actual fitted value i s 
larger than the lifetime estimated bv lKieldsen et ail |2005), 
the error bars do overlap. Estimates of the amplitude were 
also obtained, however they were rather poorly constrained, 
especially for the weaker modes. 

Simulated time series made to mimic the WIRE data 
were created in order to test the accuracy and precision of 
the fitting methods using a Monte Carlo approach. We found 
that for the most part, the fitted parameter estimates aver- 
aged over a number of realizations agreed with the input val- 
ues used to create the data. However, we did find that fits to 
the power spectrum tended to underestimate the linewidths 
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and overestimate the splittings. The bias on both of these 
parameter estimates were reduced when fitting the ACF. 

There still may be opportunities to refine this work fur- 
ther. For example, in this analysis, a fairly basic approach 
to dealing with the window function was employed, simply 
allowing for the subsequent sidebands in the models. A more 
sophisticated approach for fitting the power spectrum would 
be to convolve the spectral window with the model and fit 
that to the data. Also, because of the Weiner-Khinchine rela- 
tion, that states that the ACF is actually the Fourier trans- 
form of the power spectrum, this technique can probably be 
applied to the ACF fitting approach as well. Additionally, 
there is now a new set of WIRE a-Cen observations that 
was taken in January 2004 and lasted for around 30 days. If 
modes can be identified and fitted from this data as well, it 
will give an excellent comparison with the 1999 time series 
and should allow for a better investigation into a possible 
activity cycle. 
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APPENDIX A: VARIATION OF MAXIMUM 
POWER SPECTRAL DENSITY 

Consider a time series of length T, of a stochastically excited 
p-mode with mean power P and a lifetime r. The resulting 
power spectrum peak can then be modelled according to a 
Lorentzian function with a linewidth, A = l/(nr) . There- 
fore, assuming the modal peak is well resolved in frequency, 
the maximum power spectral density per bin, Ft, is given 
by: 

2P 
hi A 

Notice that A is multiplied by T in order to give the width in 
terms of the number of bins rather than Hz. This expression 
can be rewritten in terms of the lifetime r. 



H = 2P 



(?) 



(A2) 



However, were T to be reduced to such an extent that 
T <C t, the underlying profile would be so narrow as to 
confine all power within a single bin. In this resulting under- 
sampled regime we have conditions that tend towards an un- 
damped sine wave, where H ~ P. Therefore, a full descrip- 
tion o f H in terms of the other parameters is ijChaplin et alJ 
2003): 



H 



2P(T/t) 
P 



for T > 2r 
for T « It. 



(A3) 



Unfortunately, this does not give an adequate description for 
H when T ~ It. In order to do this we modify equation I A2I 
slightly to give: 
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Figure Al. Maximum power spectral density as a fraction of 
the full modal power determined from the maximum value of 
1000 co-added spectra (crosses). We also plot the expected values 
as determined from the oversampled (dotted line) ; undersampled 
(dashed line) and modified (solid line) expressions given in equa- 
tions EH and EH 



which can be seen to work in both the over and undersam- 
pled regimes. 

In order to test equation lA4l we produced 1000 artificial 
time series of a single p-mode signal with linewidth of 1 /iHz, 
ranging in length from 100 days down to 0.1 days. The power 
spectrum of each series was then taken and the independent 
spectra co-added to produce a smooth peak allowing the 
maximum power spectral density to be easily measured. In 
Fig. IA1I we compare these values against those predicted 
from equations IA3I and IA4I The figure shows clearly how 
the expressions given in equation IA3I diverge away from the 
true maximum power spectral density in the region T ~ 2r, 
whereas equation lA4l matches reasonably closely throughout 
the entire range. 

If one wishes to determine the power from the fitted 
linewidth, A and the maximum power spectral density then 
equation IA4I can be rearranged to give: 



H 



(f 



TA + 1 



(A5) 



(T/r) + 2' 



(A4) 



which is the equation used in Section f3. II We add that even 
though the results presented in this paper are in the range 
T > 2t, we still choose to use this more accurate expression. 
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